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Abstract 

We present a new class of thermal ratchets operating under the action of a symme- 
try breaking non-Hermitian perturbation which rectifies thermal fluctuations, and 
driven by a unbiased periodic force. The peculiar non-Hermitian dynamics which 
follows causes energy transduction from the force to the system in such a way that 
an average 'uphill' particle current is induced. We discuss physical realizations in 
assemblies of orientable particles, in itinerant oscillator models, and in problems of 
diffusion in disordered media 
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1 Introduction 



In the last few years it has been shown the existence a wide variety of trans- 
port processes at the mesoscopic level in which thermal noise plays a decisive 
role. To understand how those processes work, several physical models and 
technological implementations have been proposed [1], [2]. The peculiar effect 
of thermal noise can be illustrated in thermal ratchets or Brownian motors, 
which have the ability of extracting work from out-of-equilibrium fluctuations 
in spatially periodic systems without spatial inversion symmetry. One way 
to do this is by the combination of a ratchet-like potential which rectifies 
thermal fluctuations, and a periodic unbiased force driving the system out of 
equilibrium. 
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Our purpose in this paper is to propose a new class of Brownian motors able to 
extract work from thermal fluctuations in nonequilibrium systems driven by 
a periodic force. Unlike the ones previously introduced, the symmetry break- 
ing is due to the presence of a vortex field responsible for the existence of a 
non-Hermitian component in the stochastic dynamics. Fluctuations are then 
rectified by a non-equilibrium source instead of a ratchet-like potential. We 
will call those devices vortex ratchets. 

The paper is organized as follows. In section 2, we discuss the stochastic dy- 
namics of these systems. We formulate the Fokker-Planck equation and com- 
pute the susceptibility. In section 3, we analyze the dissipation of energy in 
the system and introduce the ratchet current. Section 4 is devoted to discuss 
some applications 



2 Stochastic dynamics 

We consider a Brownian degree of freedom, parameterized by a coordinate x, 
interacting with a thermal bath which is maintained out of equilibrium by 
the persisting action of an external drift v(x), [3]. This drift could represent, 
for example, a constant or a quenched velocity field or an external field. The 
stochastic dynamics is governed by the probability density ^(x, t) satisfying 
the conservation law 



<9 t tf(x, t) + V x • (v(x)tf (x, t)) = -V x • J^x, t) , (1) 
in which the probability current is given by 

J^x, t) = -L>V x *(x, t) + 6F(x, t)*(x, t) . (2) 

The dynamics of the Brownian degree of freedom is then governed by the 
Fokker-Planck equation 

dtV(x, t) = -V x • (v(x)W - Z)V x f) - V x • 6F(x, t)*(x, t) , (3) 

where b is the mobility, D = ksTb the corresponding diffusion coefficient, and 
F(x,t) = F e (x)A(t) a periodic force, with A(t) = \ Q e iuJt . 

In the linear response regime the formal solution of the Fokker-Planck equation 
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reads 



*(x, t) = * (x) + f \{t')^-^ c °C^ {^)dt' = * (x) + Atf(x, t) , (4) 

where C D = — V x • v + .DV X is the unperturbed Fokker-Planck operator, and 
L\ = — 6V X -F the perturbation. Moreover, \l/ (x) corresponds to the station- 
ary solution of eq. (3), [4]. 

The presence of the perturbation causes deviation in the coordinate, Ax. To 
compute that quantity, we will expand the term £i^ (x) in series of the 
eigenfunctions of the operator £ , <p n with eigenvalues a n ; n — 0, 1, .... 



A* (X) = £{ Cri 0„(x) + 4C(X)} , (5) 
n=0 

where c n are the corresponding coefficients. We obtain 

Ax(f) = J xA$(x,t)rfx = J* drxit - t)\{t\ (6) 
which defines the susceptibility xify- 

We will assume the existence of a dominant time scale governing the relax- 
ation process, corresponding to the n = 1 mode in the expansion (5). Since 
the remaining modes decay faster, we can truncate the series retaining only 
the second term. Thus, considering only contributions of the first mode the 
susceptibility is given by = Ae ait + c.c, with A defined as 

A = d J x0i(x)dx (7) 
Assuming now t Q — > — oo, eq. (6) can be rewritten as 

Ax(0 = xMA(t) , (8) 
where x( u ) is the Fourier transform of x(t) given by 



A 1 A* 1 

X ^> ~ + + ~h /3-i(a- 1) ' (9) 
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with * standing for complex conjugate. Due to the non-Hermitian nature of 
the operator, the first eigenmode is complex: a\ = R± + Hi with Ri and I\ 
being its real and imaginary parts, respectively. The remaining parameters in 
eq. (9) are (5 = Ri/h and the normalized frequency a = oj/I\. 
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Fig. 1. Non-dimensional modulus of the susceptibility as a function of a, for different 
values of the parameter /3, the smaller the value of (5 the sharper the curve. The 
resonance fades away practically for (3 ~ 10. 

In Fig. 1, we show that during the relaxation process of non-equilibrium fluc- 
tuations the susceptibility undergoes a resonant behavior when the frequency 
of the force matches the imaginary part of the first eigenvalue of the non- 
perturbed operator C Q . This behavior reveals the existence of a resonant cou- 
pling between the periodic force and the non-equilibrium source, responsible 
for the non-Hermitian nature of £ D .The implications of that coupling in the 
energy transduccion of the system will be analyzed in the next section. 



3 Ratchet effect 

Systems governed by the non-Hermitian dynamics discussed in the previous 
section may transduce the energy supplied by an unbiased periodic force into 
kinetic energy, thus inducing a net particle current. To analyze this pecu- 
liar behaviour, we will first calculate the power dissipated by the system. To 
that purpose we will apply the scheme of mesoscopic non-equilibrium thermo- 
dynamics [5]. For a system described by the probability density \l/(x, t), the 
variation of the entropy due to changes in configurations in x-space is given 
by 

5S = ~ J /i^rfx, (10) 

where \i = k B T\n^/ + U is the chemical potential, U (x, t) the potential, and T 
the temperature. The rate of change of the entropy can be obtained by taking 
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the time derivative of eq. (10), and using eq. (1). One achieves 
dS r r 

— + J Vv • V(fi/T)dx = - y -V V(/i/T)rfx , (11) 

The right hand side of eq. (11) constitutes the irreversible part of the rate of 
change of the entropy or entropy production. Consequently, the power supplied 
by the external force and dissipated into the system is obtained from eq. (11) 
by using the expression of the chemical potential. One obtains 

P F = |j^-F(x,t)rfx = F -|^(x)-(v(x))|A(t), (12) 

which defines the particle current (x) = d(x)/dt — (v(x)). To obtain eq. (12), 
we have assumed an homogeneous force and used eqs. (1) and (2). Thus, Pp 
can be interpreted as the projection of the particle current along the direction 
of the oscillating force. The quantity of interest in experiments is the time- 
averaged dissipated power, defined as 

P{u) = — / P F dt , (13) 
2ir Jo 

This quantity is not only a function of the imaginary part of the susceptibility 
as occurs in Hermitian systems. The presence of the external drift introduces 
a more complicated dependence on the moments of x. Moreover, it does not 
have a definite sign, and in general can be a non-monotonous function of the 
frequency expressing the resonant character of the energy dissipation. Since 
in general the external drift v(x) introduces a characteristic frequency in the 
system playing the same role of a vorticity, this systems behaves as a vortex 
ratchet. 



4 Applications 

Our purpose in this section is to present different manifestations of the vortex 
ratchet, as well as to indicate potential applications in different fields. 

4-1 Orientable particles 

An orientable particle of mesoscopic size in a vortex field under the influence 
of a periodic force exhibits the phenomenology discussed previously. In ab- 
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sence of the external force the orientation of the particle is characterized by 
a director vector undergoing Brownian motion which on average rotates at 
the velocity imposed by the local vorticity, u> . The most interesting situation 
occurs when the force is perpendicular to the vortex field. In that case, if 
uj > \/2uj the particle acquires an excess of angular velocity as a consequence 
of the torque it feels. In such a situation, the energy supplied by the periodic 
force is transduced into rotational kinetic energy, thus inducing a net particle 
current. This power is given by Pp — t • (ftp — |cl» ), where ftp is the average 
angular velocity of the particle, and r the torque acting on it. In this range of 
frequencies the system behaves as a Brownian motor. 

The director could represent a dipole moment oriented by a field. Suspensions 
of such dipoles in a liquid phase exhibit peculiar collective behaviours. This is 
the case of electro- and magneto- rheological fluids, ferrofluids [6], [7], [8] and 
dilute solution of rod-like polymers [9] . 

4-2 Itinerant oscillator models 

The itinerant oscillator model essentially consists of a Brownian particle with 
an orientable core coupled via an interaction potential to the shell. It was 
proposed to explain microwave dielectric absorption of polar fluids [10]. A par- 
ticular realization of the model is an inhomogeneous body under the influence 
of a vortex field and a constant force perpendicular to it. The inhomogeneity 
induces a time-dependent dipole moment. The associated torque is m(t)g x r, 
with m(t) being the dipole moment strength, g the constant external force, 
and r the orientation vector. When the dipole moment varies periodically in 
time, as in the case of the orientable particle, the system behaves as a vortex 
ratchet. Unlike the previous case, variations in the exerted torque are due to 
internal reorganizations and not to variations of the external field. The itiner- 
ant oscillator model may mimic a living cell with an inhomogeneous density 
distribution in an external field, a particle in a cage formed by other particles, 
a rod-like polymer moving in a tube [9], or a monodomain magnetic parti- 
cle whose magnetic moment undergoes fluctuations. Those systems would be 
good candidates to act as vortex ratchets. 

4-3 Diffusion through random and structured media 

A particle diffusing in a randon media advected by a steady mean-flow velocity 
v (x) in the presence of a periodic force also manifests the ratchet effect 
previously discussed. The dynamics of the particle follows from the Fokker- 
Planck equation (3 ), where the drift v(x) is now a random velocity distributed 
around the mean-flow velocity according to a Gaussian probability distribution 
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with variance T. The Fourier transform of the random drift v(k) displays both 
longitudinal and transversal correlations 

(viMvjQd)) = 2(27r) 3 r^-5(k + k') . (14) 

The transversal correlations play the role of a vorticity, introducing the non- 
Hermiticity in the dynamics. 

Assuming that the external force is homogeneous, the dissipation is obtained 
through eqs. (12) and (13), after averaging over the disorder. In Fig. 2, we have 
represented the normalized dissipated power P as a function of the normalized 
frequency a, for two levels of disorder. These levels are characterized by the 
parameter s = a 2 (Dk 2 ) 2 where D is the diffusion coefficient and a 2 = il/2k 2 T, 
with Q the volume of the system. The disorder then increases when decreasing 
s. 
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Fig. 2. Contribution of the external force to the dissipated energy as a function of 
a. Solid line corresponds to s=2, dashed line represents the same quantity for s=l. 

The figure shows how the contribution of the periodic force to the total dis- 
sipated power exhibits a minimum when the frequency matches the charac- 
teristic frequency of the system, thus revealing the resonant character of the 
dissipation. Due to the characteristics of this system, when the power is neg- 
ative the particle current is positive. The figure also shows that the power or, 
in view of eq. (12), the current is positive for frequencies around the resonance 
frequency, which manifests that the velocity of the particles is larger than the 
average drift. Therefore, in those conditions the system acts as a Brownian 
motor. 

A similar phenomenon occurs in a spatially periodic two-dimensional pattern 
of triangular vortexes when a two-dimensional oscillating force is applied, at 
sufficiently high Reynolds number. In this conditions, a large-scale current 
appears. This effect is accompanied by a reduction of the dissipation in the 
system due to the induction of a negative eddy viscosity[ll]. 
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